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CELFE/ NASTRAN CODE FOR THE ANALYSIS OF STRUCTURES 


SUBJECTED TO HIGH VELOCITY IMPACT 
by C. C. Cham is 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

ABSTRACT 

The CELFE (Coupled Eulerian Lagrangian Finite Element)/NASTRAN 
Code three-dimensional finite element code has the capability for analyzing 
of structures subjected to high velocity impact. The local response is pre- 
dicted by CELFE and, for large problems, the far-field impact response is 
predicted by NASTRAN. The coupling of the CELFE code with N AST RAN 
(CELFE/NASTRAN code) and the application of the code to selected three- 
dimensional high velocity impact problems arc described. 


INTRODUCTION 

A three-dimensional finite element computer code was recently developed 
for determining the response of components and structures which are subjected 
to high velocity impact. The schematic in figure 1 depicts a general class of 
high velocity impact problems for which a three-dimensional finite element 
code is required. This code is based on the coupled Eulerian- Lagrangian for- 
mulations (mode) and has been given the acronym CELFE (Coupled Eulerian 
Lagrangian Finite Element). CELFE is structured to handle high velocity im- 
pact structural response in the presence of geometric and material nonlineari- 
tics, material flow, and anisotropic material behavior. CELFE has been cou- 
pled with NASTRAN in order to: (1) take advantage of already available capa- 
bilities in NASTRAN to solve large systems of equations which require out-of- 
core storage, and (2) to determine the f r " ?!d structural response of large 

structures subjected to local high velocity impact. 

The code was developed by the Lockheed Missiles and Space Company, 
Huntsville Research and Engineering Center, Huntsville, Alabama, under con- 
tract NAS3-18908 to NASA Lewis Research Center (LeRC). It is completely 
documented in references 1 to 3. The code was installed on the LeRC com- 
puter (UNI VAC 1110) in late January 1978. Presently, our effort has been 
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devoted mainly to becoming familiar with the code capabilities and to correcting 
system-associated errors. It is anticipated that the code will be available for 
public distribution through COSMIC (University of Georgia, Athens, Georgia 
30601) by early 1979. 

The objective of the present paper is to provide a brief description of the 
CELFE/NASTKAN Code and it. 1 capabilities, and discuss selected problems 
which have been solved using the coupled computer program. Detailed descrip- 
tions of all aspects, that is, formulation, programming and user's manual, are 
given in references 1 to 3. The notation used in the disc 1 •''•ion is defined when 
it first appears and it is also summarized in the appendix for convenience. 

COUPLED El' LEIUAN- LAGRANGIAN FINITE 
ELEMENT (CELFE) CODE 

The CELFE code used in the CELF.3/NASTRAN program has the following 
general capabilities: 

First, it simulates the transient state during the impact process and right 
after the process is completed; second, it computes the dynamic response of 
the entire structure for specified times after the impact process is completed; 
and third, it predicts various phenomena which occur during the impact process 
such as rebounding, sliding, and penetration. Also, CELFE traces the move- 
ment of the free surface, and the failure front. The specific capabilities of 
CELFE are summarized in table I. 

The CELFE code is divided into three pails: the first pail contains a 
finite element module to simulate the dynamic behavior in the impact Eulerian 
zone; the second pail consists of a regular finite . . .* procedure to analyze 

the dynamic response in the near field (Lagrangian) zone; and the third part 
consists of a coupling (interfacing) procedure to integrate the above two sepa- 
rate pails to form a complete system, that is, the CELFE system. The 
CELFE system has two options for computational purposes. The in-core 
CELFE code option (for small problems) and the coupled CELFE/NASTKAN 
rc ’o option (for large problems). The in-core option is a subset of the cou- 
pled option to allow more flexibility in utilizing the computer program. 


STRUCTURAL MODELING AND FORMULATION 

The CELFE Code treats the entire structure as partitioned into an impact 
zone and a structural (lagrangian) zone. The impact zone co.irists of an 
Eulerian zone (E zone) and a transition zone (E-L zone). The Lagrangian 
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zone (L zone) can be divided into an in-core and a N AST RAN zone according to 
the opt lonn of the In-core CELFE code and the coupled CELFE/NASTRAN code. 
The zones are denoted L t . and L^, respectively, as depicted schematically in 
figure 2. Since the in-core option is a subset ol the coupled option, the La- 
grangian zone (L zone) may be composed solely ol L, when the in-core option 
is used; or composed solely of L v or mixed L . and L v when the coupled 
option is used. This allows the user to arrange his problem to optimize com- 
putational efficiency. 

In the entire impact zone, together with the L c . zore, three-dimensional 
solid elements are used, in the remaining portion of the structure, that is, 
the L^- zone, any appropriate elements (depending on the nature of the struc- 
ture) can be applied. In the impact zone (K zone and E-L zone), the dynamic 
behavior is handled by the three-dimensional impact element using a hydro- 
elasto- viscoplastic model with a general moving coordinate system. The pri- 
mary variables in this zone consist of the density, p, momentum, pV jf j 1, 

2, 2, total energy , pc, and deviatorie stresses, S 1)t i, j = 1, 2, 3. The remain- 
ing variables (pressure P, stresses o- jt and displacements u.) belong to the 
secondary set. In the Lagrangian zone (L^. and L^. zones) conventional struc- 
tural dynamic analyses are applied. In this study, the displacement method 
was adopted, and thus the displacements (both the lateral and rotational dis- 
placements) were chosen to be the primary variables. The governing equations, 
coordinate systems, and types of elements in the various zones are different. 

To accommodate these, the coupling procedure provides two mechanisms; 
namely, an in-core procedure to couple the E-L and L . zones (or L». zone in 
case of empty L c ); and an interfacing procedure to couple the main CELFE 
program with NASTRAN. The coupling variables for the present procedure 
contain lateral displacement and velocity fields. 

The governing equations of the finite element formulation are summarized 
in table 11. The associated variables, types of elements, and related aspects 
concerning the above discussions for various zones are summarized in table 111. 


SOLUTION PROCEDURE USED BY THE COMPLETE CELFE SYSTEM 

The CELFE system employs a generalized two-step time integration 
scheme, integrated with respect to time. For the coupled CELFE/NASTRAN 
option, it may not be necessary to include the NASTRAN part. In the present 
code, it is included for the sake of simplicity in coding the interface procedure. 
The solution procedure includes both the in-core CELFE option and un coupled 
CELFE/NASTRAN option. The solution procedure is summarized in the flow 
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chart shown in figure 3. The key computations in this flow chart are as follows: 

1. Read in parameters defining the problem and start integration with re- 
spect to time (boxes 1, 2, and 3). 

2. For each time step: 

a. Predict the solutions of primary variables in various zones using 
the results obtained in the previous time step (box 5). 

b. Construct the global system matrix for the entire structure in 
terms of coupling variables, and solve the resulting equation using the 
solution obtained in the previous time step (box 8). 

c. I pilule the secondary variables in various zones using the above 
predicted values of the primary variables (box 9) and repeat steps a and b. 

d. I pdate the mesh in the E-L zone and/or L c zone (box U). 

e. Correct the results by repeating steps a to d using the above pre- 
dicted solutions together with the results obtained in the previous time 
step. 

1. Test the corrected results tor the criteria governing projectile 
rebounding, sliding, and/or penetration (box 10). 

3. U the test shows the penetration continues, repeat the steps stated in 
item 2 above for the next time step. Otherwise, 

a. II the projectile rebounds from the target, and/or if disturbances 
in the structure are damped to a certain user selected fraction of the 
original values, the program regards the entire structure as a Lagrangian 
zone and the dynamic response is handled accordingly (box 7). 

b. If the projectile slides over the target due to oblique impact, the 
program considers the mesh for the entire structure to be Lagrangian 
(box 7). The steps described in item 2 above must be repeated. 

In order to maintain the accuracy of computations in the vicinity of the 
impact point, as well as to have a manageable in-core storage, the dimensions 
of the impact zones in this code are appropriately chosen a priori. At t = 0, 
the entire impact zones are identical with the E-L zone, where the zone is ini- 
tially discretized into Lagrangian mesh (fig. 2). For ' * U, when failure is 
detected at certain nodes, the mesh containing these nodes is turned to Eulerian. 
Hence, the impact zone, whose dimension has to be specified a priori, must be 
sufficiently large to contain the failure region up to the impact process comple- 
tion. If excessively large impact zones are assigned, on the other hand, it will 
require unnecessary computer storage and computational time. 
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PRCCEDl HE FOR COl PUNG CELFE WITH N AST KAN 
(CELFE/NASTKAN) 

The coupled CE LFE/NASTRAN code is organized to solve large structural 
systems as shown in the How chart in figure 4 in order to handle either or both 
of the to How mg: 

1. Utilize the matrix equation solver SOLVE of N AST KAN to obtain the 
solution in the impact zone; and/or 

2. Couple the Eulerian and Lagrangian modes for the entire structural 
system. 

I tilization o' N AST RAN Subroutine SOLVE 

Symbolically, the governing equation (first equation, table li) can lx* re- 
written for each element in the impact zone as 

Av> = B (1) 

where «.■ is a vector of the 14 primary variables (density, internal energy, 
three velocities, and nine stresses); A represents geometry, material, inertia 
and dynamic relationships (Eulerian formulation) between the primary variables 
at the current time step; and B contains A and c from the previous time 
step (inertial forces) and the associated dynamic forces. 

Equation (1) can be solved using the N AST RAN SOLVE module for each 
primal*} variable. In order to optimize the time required for the interfacing 
procedure, the following scheme is used: 

1. At each time step, individually compute A and B for each element, 
and for each primary variable. These element matrices are stored on a disk 
file in INPUTT2 matrix format (boxes 0, 7, and 9). 

2. Store all relevant data, including the solutions of the previous time step, 
and those of the secondary variables, etc. into disk files (D.,, D. { , and P^) 
before calling EXIT from CELFE. 

:L Using N AST KAN, read (box 11) all the data of element matrices A and 
B (box 12), and sum the element matrices for each primary variable (box 14). 
Then solve equation (1) using SOLVE (box 14). 

4. Write the solutions for ^ on a disk file using OUTPUT2 matrix format 

0 > 4 >* 

5. I sing CELFE, read the output disk and make further computations 
(box 5). 
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Coupling CELFE with NASTRAN (CELFE/NASTRAN) 


The governing equations for this part are: 

1. CELFE (impact) zone tor each element (denoted by superscript) 




( 2 ) 


where K is the equivalent stiffness matrix, u the displacement vector, and 
G the equivalent forces all for the current time step. 

2. Lagrangian and/or NASTRAN zones 

|[p 2 M| L+ [pB] L+ |K| L j <L, {6} ( *> = {G}£> (3) 

Summing equations (2) and (3) over e = 1, 2, . . . , N, yields 

|[p 2 .MJ + [pBJ + [Kjj {6} ={G} (4) 

N is the total number of elements in the entire structure (including the projec- 
tile), p denotes time derivative (d/dt), and 6 denotes global displacements. 

The actual coupling variables consist of lateral displacements and veloci- 
ties, Uj and Vj, j = 1, 2, 3, in the impact zone. The lateral displacement 
field consists of six degrees of freedom (three lateral and three angular dis- 
placements). 

The procedure for the solution of this part is similar to the in-core version 
with some additional steps to couple CELFE with the NASTRAN as follows: 

1. Compute all element matrices (eq. (2)) at each time step and store them 
in mass storage formated for the NASTRAN module INPUTT2 (box 7 and D.,, 
flow chart in fig. 4). 

2. Compute the load vector G and store it in disk file (box 9; P^). Store 
all relevant data in disk files (L^ and D.j) before calling EXIT from CELFE 
(box 10). 

3. Call NASTRAN (box li) 

a. Read all element stiffness matrices from the disk file and add the 
matrices into NASTRAN structural stiffness matrix K in equation (3) to 
form the coupled system matrix K (box 13). 

b. Add the inertia and damping terms, together with the load vector 
(box 12) composed of those read from the breakpointed disk file for CELFE 
grids and form equation (4) (box 13). Solve the system and store the results 
in OUTPUT2 format on a disk file (box 13). 

c. Return to CELFE for further computations (boxes 15 and 4). 
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The procedures utilizing the equation solver module SOLVE in NASTRAN 
lor the impact zone, and the coupling of CELFE with NASTRAN ran be used 
either separately or simultaneously. The choice depends on the specilic 
problem under consideration. 


ANALYSIS OF A THREE-DIMENSIONAL IMPACT PROBLEM 

The results of a three-dimensional impact problem are discussed in this 
section to illustrate the use of the CELFE/NASTRAN program. The given 
conditions of the example are as follows: 

Target: boron/epoxy (AVCO 5505) plate with dimensions 9.75 by 8.75 by 
0.U83 inch 

Projectile: silastic with dimensions (0.30)^ in 

Impact velocity: 030 ft/sec normal to the plate 
The geometric configuration is depicted in figure 5, and material properties 
are listed in table IV. The boron/epoxy target was unidirectional with fibers 
parallel to the y-axis. 

Due to the symmetry, only one quadron of the plate needed to be consid- 
ered. The finite element mesh of the model is shown in figure 0. The impact 
zone consisted of the first 33 elements (nodes 1 to 70). Elements 30 to 33 
(nodes 02 to 70) constituted the transitions zone (E-L zone). 

In the Lagrangian zone all elements are NASTRAN general quadralateral 
elements (L v zone). The L (i zone was assigned to be empty for simplicity. 

The mesh is illustrated in figure 7 . To ensure displacement compatibility 
across the boundary, multipoint constraints (MPC) were employed to connect 
the CELFE and NASTRAN zones. A NASTRAN grid point was added to each 
CELFE midpoint grid lying on the E-L zone. This E-L zone has now become 
the interface of the CELFE and NASTRAN zones due to the empty L Q zone. 

The NASTRAN grid points added were 163, 1(36, 169, 172, and 175, respec- 
tively, to the CELFE grids 63, 66, 69, 72, and 75 (fig. 7). The midpoint 
grids were coupled by MPC's for the tl ree translational degrees of freedom. 
Rotational displacements were transmitted from the quadralateral plate ele- 
ments to the CELFE isoparametric elements by using rigid bars connected 
between the midpoint node (node 163, fig. 8) and the top and bottom nodes of 
the solid (nodes 64 and 62, respectively). 

Results for the example problem are illustrated in figures 9 to 11. In fig- 
ure 9 the pressures developed in the plate midplane nodes (5, 8, and 11, fig. 6) 
are plotted as a function of time. The pressures build up rapidly and are ap- 
proximately the same for all three nodes. The three normal stresses (a , 

AA 

a, , and u_ ) are plotted along an axis passing through the impact center after 
yy zz 



0. 07 microsecond (fig. 10). All the stresses are compressive, increase rapid- 
ly In the contact region, and decrease very rapidly away from this region. In 
figure 11 the failure propagation is shown as a function ol time. Failure occurs 
almost at the initial stages oi the impact event. The results obtained from this 
example problem demonstrate how the CELFE/NASTRAN code can be used to 
describe the local and far field response of structural parts subjected to high- 
velocity impact. 

The computer CPU times (t'MVAC lies) lor this problem were 0 minutes 
pe r time step in the impact zone and do minutes per time step for CELFE/ 
NASTKAN. These times are relatively large and speed improvements are re- 
quired to make the code practical lor parametric and sensitivity studies. 


SUMMARY 

The important characteristics oi the CELFE/NASTRAN code are summa- 
rized below: 

1. A three-dimensional computerized cajiability has been developed (CELFE/ 
NASTKAN) which can be used to describe the local and far field responses of 
structures subjected to high velocity impact. 

2. The CELFE part of the capability accounts lor local geometry and mate- 
rial nonlinearities, anisotropic and viscoplastic material behaviors, material 
flow and penetration, normal, oblique and sliding impact. 

d. The coupled CELFE/NASTRAN code can handle far field material and 
geont et ry non 1 i nea r it ie s . 
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APPENDIX - SYMBOLS 


A array of coefficients (geometric, material, inertial, flow) of inter- 

relationships between the primary variables in the Eulerian mode 

(eq- U)) 

A array of displacement shape-function relationships 

ID 

B vector of inertial and dynamic forces in the Eulerian mode (eq. (1)) 

B fp array of relationships between the derivatives of the displacement 

shape functions and those of the primary variables 

CELFE Coupled Eulerian Lagrangian Finite Element 

C f vector of relationships between velocity and displacement shape 

functions 

c subscript - compression 

D mass storage disk files used in code for storing arrays 

E Eulerian zone 

E-L Eulerian- Lagrangian coupling zone 

Ej modulus of elasticity, subscript? •* * Jte direction 

e internal energy 

G equivalent element forces for the current time step (eq. (2)) 

shear modulus, subscripts denote direction 
1NPUTT2 NASTRAN module for input data 

i index 

index 

K equivalent element stiffness matrix for the current time step 

(eq- (2)) 

L lagrangian zone 

L e CELFE zone 

NASTRAN zone 

M equivalent element mass matrix (eqs. (3) and (4)) 

NASTRAN NASA STRUCTURAL ANALYSIS 

n • script - current iteration cycle 
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OUTPUT2 


P 


P 

S 



s 


( 


T 

t 

t.U 

V, v 


x.y.z 


1 , 2,3 


3-D 

6 

'l 

P 

pe 

pv 

% 


N AST HAN module lor NASTHAN output data 
pressure; mass storage disk files for load vectors 
derivative with respect to time (d/dt) 
subscript, shear 
deviatoric stress tensor 

.. . . composite strength (subscripts denote direction and 

sense) 

subscript - tension 
time 

time step 

displacement 

velocity 

orthogonal coordinate reference system; subscripts denote coordi- 
nate directions 

subscripts where 1 is along the fiber direction, 2 transverse, 
and 3 through the thickness 

three-dimensional 

global displacement 

Poisson's ratio, subscripts denote direction 

density 

total energy 

momentum 

stress field 

vector of primary variables in the Eulerian mode (eq. (1)) 
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TABLE I. - SPECIFIC CAPABILITIES OF THE CELFE CODE 

1. Predict normal and oblique impact 

Penetration process 
Sliding process 
Rebound process 

2. Penetration process 

Dynamic analysis during penetration process 
D>namic response ol the structure alter penetration pro- 
cess is completed 

3. Relevant criteria in CELFE system 

Failure criterion for composites 

Plastic and viscoplastic yielding for isotropic materials 
Sliding criterion - modified coulomb's law 
Rebound criterion 
Degeneracy due to damping 

4. Large structure impact via CELFE/NASTRAN coupling 

TABLE 11. - GOVERNING EOL ATIONS OF THE 
FINITE ELEMENT FORM ELATION 
[See appendix for symbol definitions. J 
Finite element analog of impact zone: 

(l A r S J * f l«w) { u js +1 ’} (t A rsl + I B rs l) {<} ♦ ^!C r } 

Finite element analog of Lagrangian zone: 

[lP 2 M] + lpI3] + [Kjj {6} ={G} 

Coupling variables: 


Uj, V-; J * 1, 2, 3 



TABLE III. - CLASSIFICATIONS OF VARIOUS ZONES IN HIGH VELOCITY IM PACT 


(See appendix for symbol definitions.] 


Classification of 
zones 

Impact zone 

Lagrangian zone 

E-zone 

E-L zone 

L c zone 

L v zone 

Eulerian zone 
(rrnpiy at t = 0) 

Transition zone 

In-core Lagrangian 
zone (may be empty) 

N AST RAN zone 
(may be empty) 

Primary- variables 


j = 1,2,3 

Secondary- variables 

u j» P,cr ij ; = 1,2,3 

Uj.Oj; j = 1,2,3 (Uj = Vj) 

Type of elements 

3-D isoparametric elements 

3-D isoparametric 
elements 

Any compatible 
elements 

Coupling variables 

Vj, Ujj j — 1, 2 t 3 
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TABLE IV. - PROPERTIES EUR A BORON AVCO 5505 PLATE 
AND A SILASTIC PROJECTILE 
[Subscripts defined in the appendix. J 


Target 


Projectile 


Material 


Boron AVCO 5505 


Silastic 


0.073 


0.0512 


Elastic modulus, psi 


E (U *= 29.2*10 u 
E l22 " E i 33 315xi0< 


1.79 lo 


Shear modu’us, psi 


Poisson’s ratio 


G u2 ^ G ( 13 =• 0.78- 10 u 0.0. 10' 

u (23 = o.a-ao 6 


"tl2 = ''(U “ 017 


^23 = 0 - 53 


Uniaxial failure stresses, ksi S^^*19!I.O 


S £11C = 232 0 
S (22T S (33T ~ 8,1 

S £22C " S (33C = 17,9 
S (12S S £ 13S " 9,1 
S ( 23S ~ 8,9 



Yield: 27027 



















IMPACTING 

OBJECT 

iPROJECTUEi 


LAYIREO 

ANISOTROPIC 

MEDIUM 


POINT Of 
IMPACT 


im 


GHE> 


)*D LOCAL 

IMPACT 

ElfMENT 

IMPACT REGION 
3-D ElfMENT 

SIRUCTURAl 
COMPONENT 
iTARGETl im.vl 



OTHER 

ELEMENTS 


Figure 1. • Impact conditions and idealizations. 
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Figure 2 . - Typical CELFE/NASTRAN finite element model. 




Figure 3. - Flo* char! of CElfE solution procedure. 














Figure 4. - Flow chart Interfacing procedure tor CELFEJNASTRAN. 
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Figure 11. • Propagation of failure front with respect to time. 



